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ABSTRACT 

The high-mass X-ray binary RX J0146.9+6121, with optical counterpart LS I+61°235 (V831 
Cas), is an intriguing system on the outskirts of the open cluster NGC 663. It contains the 
slowest X-ray pulsar known with a pulse period of around 1400s and, primarily from the study 
of variation in the emission line profile of Ha, it is known to have a Be decretion disk with 
a one-armed density wave period of approximately 1240d. Here we present the results of an 
extensive photometric campaign, supplemented with optical spectroscopy, aimed at measuring 
short time-scale periodicities. We find three significant periodicities in the photometric data 
at, in order of statistical significance, 0.34d, 0.67d and O.lOd. We give arguments to support 
the interpretation that the 0.34d and O.lOd periods could be due to stellar oscillations of the 
B type primary star and that the 0.67d period is the spin period of the Be star with a spin 
axis inclination of 23^g° degrees. We measured a systemic velocity of —37.0 ± 4.3 km s _1 
confirming that LS I+61°235 has a high probability of membership in the young cluster NGC 
663 from which the system's age can be estimated as 20-25 Myr. From archival RXTE ASM 
data we further find "super" X-ray outbursts roughly every 450d. If these super outbursts are 
caused by the alignment of the compact star with the one-armed decretion disk enhancement, 
then the orbital period is approximately 330d. 

Key words: accretion: accretion discs - stars: high-mass X-ray binaries - stars: neutron stars 
- stars: close binaries - stars: Be stars 



1 INTRODUCTION 

High-Mass X-ray Binaries (HMXBs) are interacting binary stars 
which emit X-rays as a result of mass transfer and accretion pro- 
cesses. Be X-ray binaries (BeXs) consist of a rapidly rotating main 
sequence Be star with an equatorial decretion disk, which gives rise 
to the optical emission lines, orbited by a compact object that is 
usually a neutron star. The other major class of HMXBs is com- 
posed of the supergiant X-ray binaries (SGXs) that consist of a 
more evolved OB star closely orbited by a neutron star (or in some 
cases a black hole, like Cyg X-l). SGXs tend to have short (a 
few days) orbital periods while the orbital periods of BeXs are 
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considerably longer. BeXs are generally transient X-ray sources 
with outbursts occurring when the compact star passes periastron 
( IOkazakil200TI) . All HMXBs are relatively young, rapidly evolving 
systems in which one component, in most known cases, has already 
evolved to a supernova and subsequently to the current compact ob- 
ject. Mass redistribution within the system together with mass loss 
may prevent the supernova, however, and leave a white dwarf as 
the compact object. The H MXB 7 Cas, for example, is suspected 
of co ntaining a white dwarf jLopes de Oliveira et al.|[2006l : lHaberll 
1 19951) . Since a significant fraction of all stars form as binary stars, 
a complete understanding of the evolution of high-mass stars re- 
quires an understanding of their evolution in the presence of mass 
transfer in the interacting binary situation. 

The physical properties of many HMXB systems, in partic- 
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ular their orbital period, are known from their X-ray behaviour. 
But the orbital periods of many Galactic systems (approximately 
half of those w ith known optical counterparts) remain unknown 
dLiu et al, I l2006h . We have recently initiated a long-term program 
to obtain time-series optical photometric data for bright HMXBs in 
an effort to measu re periodicities and determine the orbital periods 
dSartv et alj|2007t) . Here we report the results of one of our first 
observing campaigns aimed at RX J0146.9+6121/LS I+6T235. 
The data were obtained by observers associated with the Ameri- 
can Association of Variable Star Observers (AAVSO) in response 
to our call to observe HMXBs on an on-going basis. In addition 
to the photometry, we also made optical spectroscopic radial veloc- 
ity measurements using the Dominion Astrophysical Observatory's 
Plaskett telescope. Significant photometric and radial velocity vari- 
ation was found which we report and interpret here. We further 
interpret the optical data alongside publicly available X-ray light 
curves from the All Sky Monitor (ASM) on board the Rossi X- 
ray Timing Explorer {RXTE) satellite and previously published in- 
formation to make an estimate of the orbital period. This paper is 
organised as follows. The optical and X-ray observations are de- 
scribed in Sect. [2] The results of period searches in the optical and 
X-ray light curves are presented in Sect. [3] Interpretation of the re- 
sults is given in Sect. [4] where we present evidence for /3 Cephei 
type pulsations in the optical component of the system, confirm its 
membership in the star cluster NGC 663 and argue that the orbital 
period is about 330d. We conclude the paper in Sect. [5] 



2 OBSERVATIONS 
2.1 Optical photometry 

Three observing seasons of optical photometric data have been ob- 
tained of LS I+61°235 covering the ends of 2005, 2006 and 2007. 
These data were obtained primarily in the Johnson V band by the 
observers listed in Table[T] where they are identified by the observer 
initials assigned to them by the AAVSO. The complete V band 
dataset is shown in Fig. [2] where it may be seen that the first two 
seasons were covered exclusively by HUZ. Cousins Ic band data 
are shown in Fig.[3]and a limited set of Johnson B—V colour obser- 
vations were obtained by LDJ. Other observers joined the HMXB 
observation program in th e third season in re sponse to a general call 
for HMXB observations dSartv et alj|2007l) and to AAVSO Alert 
Notices 348 and 354. All the photometric observations are avail- 
able from the AAVSO International Database. The total numbers of 
observations are shown in Table [2] The observations were initially 
retrieved from the AAVSO International Database and then the ob- 
servers were individually contacted to resolve discrepancies and er- 
rors, if any. This "validation" step is necessary to restore some ho- 
mogeneity to an otherwise relatively inhomogeneous dataset. All 
observers reported an estimate of the signal to noise ratio (S/N) 
of each measurement, usually based on a formal estimate of the 
Poisson shot noise. All of the photometry reported was differen- 
tial photometry with the magnitudes being determined relative to 
the comparison stars labelled in Fig.Q] A variety of methods were 
used by the observers to reduce the data using aperture photom- 
etry and software like MaxImDL, Aip4Win or MiraAP. None of 
the data were transformed to the standard Johnson/Cousins photo- 
metric system and therefore are instrumental magnitudes. Since our 
main aim is to determine periodicities from the light curve, instru- 
mental magnitudes are as valuable as properly standardised data 
provided there are no large unaccounted offsets between the sub- 
sets. The majority of our data were obtained independently by two 
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Figure 1. Comparison stars used by the AAVSO observers for the dif- 
ferential photometry reported here. The target star, RX J0146.9+6121/LS 
I+61°235, is marked by the two dashes. The data are based on measure- 
ments made at the Sonoita Research Observatory (SRO) by Arne Henden 
and the constancy of the stars has been verified from data obtained by ob- 
servers HUZ and CTE. This image is based on the POSS2/UKSTU Red 
images from the STScI Digitized Sky Survey. The field of view is 0.25°. 



Table 1. List of observers and their equipment. 



Observer 


Telescope 


CCD 


HUZ, R Huziak 


30-cm SCT 


SBIG ST-9XE 




30-cm SCT 


STL- 130 IE 


CTE, L Catalan 


25-cm Newtonian 


KAI-2020M 


LDJ, DJ Lane 


28-cm SCT 


SBIG ST-9XE 


CTX, TR Crawford 


30-cm SCT 


SBIG ST-9XE 


PXR, RD Pickard 


30-cm SCT 


Starlight Xpress SXV-H9 


CPE, P Closas 


15-cm refractor 


Starlight HX5 16 


OSC, SC Orlando 


20-cm SCT 


SBIG ST-8XE 


GTZ, TA Grzybowski 


25-cm SCT 


SBIG ST-402 



observers. The photometric values reported by those two observers, 
at times of simultaneous observation, show excellent agreement. 



2.2 Optical spectroscopy 

Optical spectroscopy at a resolution of approximately 0.4 A (10 
A /mm, A/AA ~ 1 1000) was obtained during two observing runs 
with the Cassegrain spectrograph on the 1.85-m Plaskett telescope 
at the Dominion Astrophysical Observatory (DAO) in Victoria, 
British Columbia, Canada. The two observing runs were from Au- 
gust 20 to September 3, 2007 and from December 3 to December 
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Figure 2. The complete V band light curve. The first season covers HJD 2453583 to 2453808, the second season covers HJD 2454019 to 2454162, and 
the third season covers HJD 2454347 to 2454457. No long term changes greater than about 0.05 magnitude are seen over the approximately 1000 days of 
observations. The two "V" symbols mark the times of the two largest X-ray outbursts during the three observing seasons as seen by the RXTE ASM and shown 
in Fig.|5la). 
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Figure 3. Cousins 1q band photometric data provided by observer CTX. 



Table 2. List of total numbers of photometric observations by the AAVSO 
observers broken down by filter. 



2.3 RXTE ASM X-ray data 

Archived data from the RXTE ASM were used to identify out- 
bursts that occurred before and during our optical observing runs. 
The RXTE ASM is composed of three Scanning Shadow Cameras 
(SSCs) that perform sets of 90s pointed observ ations (dwells) cov - 
ering about 80% of the sky every ~90 minutes jLevineetalJl996t) . 
Definitive ASM data for RX J0146.9+6121 were downloaded from 
the RXTE Guest Observer Facility. The data used were from the 
dwell by dwell compilation. Each raw data point represents the 
fitted source flux from one 90 second dwell. Data from all three 
SSCs were used and represent nominal 2-10 keV rates in ASM 
counts/second (c/s). Nominally, the Crab nebula flux is about 75 
ASM c/s (when the source is at the centre of an SSC field of view 
and all 8 anodes are operational). We retained only the "three- 
sigma" detections based on the reported variance for each data 
point. The ASM data are shown in Fig. [5] 
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12, 2007. The spectrograph was nominally set to receive the band 
from 4300 to 4550 A with a 0.005-in slit in the "21181" configura- 
tion, meaning that a 21 -in focal length camera, 18-hundred grooves 
per mm grating and the 1st spectral order were used. The CCD 
was a SITe-2, with 1752x532 15 ^im pixels. Exposure times var- 
ied from 900 to 2400s. The spectra were extracted and wavelength 
calibrated against the emission lines from an Fe-Ar lamp using the 
IRAF do si it task. 

A total of 1 1 spectra were obtained as shown in Fig. [4] Of 
those spectra, 3 were quite noisy because of passing clouds. Radial 
velocities were measured from these spectra as described in Sect.[3] 



3 RESULTS 

3.1 X-ray outbursts 

The previously repor ted outburst of 1997 dMereghetti et alj|200fj : 
lReig & Roche] 1 19991) can be seen in the X-ray light curve at 
HJD 2450769. Subsequent outbursts are easily identified, with the 
largest outburst occurring at HJD 2454083. The large outburst at 
HJD 2454083 (see Fig. [5} corresponds to a possible brightening of 
the V band light curve in the second observing season peaking at 
approximately HJD 2454050 (2006, see Fig.O and appears to be 
the peak of X-ray activity lasting about 40 days. The peak of the 
2006 V band brightening appears to precede, by a few tens of days, 
the X-ray peak. The peak in V band brightness seen at about HJD 
2453630 in the first observing season (2005) precedes, by roughly 
80 days, the X-ray outburst marked by X-ray flux of about 7.5 ASM 
c/s (~0.1 Crab) at HJD 2453710 and 2453770. In this case the V 
brightness appears to have faded considerably (by about 0.03 mag) 
by the time the X-ray outburst occurs. 

A formal se arch for periods in the ASM data was done using 
the period04 i Lenz & Bregen 120051) so ftware, a Lomb-Scargle 
periodogram jLomblll976l : IScarglelll982T) and a phase dispersion 
minimisation (pdm) technique i Stellingwer i 1 9781) . All methods re- 
port a strong 184d period although the third subharmonic of that 
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I+61°235 obtained from two DAO observing runs. The spectra are dominated by four lines: H-y (A4340), Sill (A4430) and He I 
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Figure 5. (a) RXTE ASM 2-10 keV X-ray light curve. The dashed line 
shows the mean background level, which is essentially subtracted from the 
three-sigma detections, and the time span of our photometric observations 
is shown. The major outbursts are separated by about 450d with smaller out- 
bursts occurring approximately every 90d. (b) The window function for the 
sampling pattern underlying the ASM light curve. There is no large concen- 
tration of spectral power in the window function at a period of 1 84d (0.0054 
cycles/day) although a small local maximum occurs there (marked with the 
"V" symbols). 



period (552d) is stronger in the pdm analysis. The window func- 
tion, representing the Fourier transform of the light curve sampling 
pattern, shows no structure (strong peaks) that could cause a 184d 
period to appear in the periodograms (Fig. [5jb)). However, folding 
the X-ray light curve at 184d reveals no obvious pattern. A visual 
inspection of the ASM light curve reveals an outburst roughly ev- 
ery 90d with every fifth outburst (every 450d) being of much larger 
amplitude. 



3.2 Ic band and B — V colour 

Referring to Fig. [3] the Ic band data shows a light curve shape very 
similar to that seen in the V band data (see below). At this point 
however, we cannot be certain that the Ic light curve varies with 
the same frequency and phase as the V light curve. Simultaneous 
(or nearly so) observations in Ic and V are required to assess the 
relationship between the two light curves or, equivalently, if there 
are changes in the V — I colour. 

Three nights of B — V data from observer LDJ show signs of 
variability, indicating that the amplitude and/or frequency and/or 
phase of the B and V light curves are different. Again more 
multi-colour photometry is required before anything positive can 
be claimed about colour variations. 



3.3 Photometric period analysis of the V band data 

Iterative Fourier analysis of the V band data, shown in detail in 
Figs.[6]and|7]for the s econd and third seaso ns, was done using the 
software period04 dLenz & Breger|2005i) in which the data were 
prewhitened with the dominant frequency between iterations. Only 
the data from observers HUZ and CTE were included in the Fourier 
analysis because these data consistently had millimag precision and 
were compatible without transformation to a standard photometric 
system. The analysis was repeated three times; once for the whole 
dataset (840 days), once for the second season (160 days), and once 
for the last part of the last season (30 days). In addition, each of the 
three analyses was done twice, once without weighting the data and 
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a second time weighting the data with weights w = 1 /a 2 where 
a is the reported error. The frequencies found from the iterative 
weighted Fourier analysis are shown in Table [3] 

Selected periodograms associated with the entire dataset (all 
three observing seasons) are shown in Fig. [8] For the entire dataset, 
both the weighted and unweighted analysis found a low frequency 
component with the highest peak with a period of almost exactly 
840 days, i.e. the time-span. This corresponds exactly to the ex- 
pected effect of the single mean brightness offset in the first half of 
the central chunk of data (second season). The 1.49 cycles/day (c/d) 
frequency is roughly half (not significantly different from half) that 
of the prominent period at 2.92 c/d and folded phase graphs at 
those two frequencies produces plots (Fig. [9]( that are very simi- 
lar to directly observed all night time series (see Figs. [6] and [7). 
The 0.18 and 1.21 c/d periods appear to be spurious because of the 
distorted shape of the folded curves, sampling window effects, a 
non-periodic mean brightness shift, etc. The whole spectrum above 
~5 c/d can be very well approximated with a 1// noise which may 
have an astrophysical origin. The 9.679 c/d peak is, however, very 
narrow, very strong, and much higher than the noise in the sur- 
rounding area of the spectrum. A typical method for assessing the 
signal to noise (S/N) ratio for a peak in the spectrum is to calculate 
the ratio between its amplitude and the mean amplitude of the sur- 
rounding peaks with a S/N > 3.5 to 4 being a commonly adopted 
limit of significance (very roughly representing, via the central limit 
theorem, a probability of arising by chance of <5%). The 9.679 c/d 
peak is significant by this definition having a S/N in the range 4.5 
to 7.6, depending on width of the frequency bins. The frequencies 
of 2.92 c/d (and/or 1.49 c/d) and 9.68 c/d are therefore judged to 
be real from the analysis of the whole dataset. There was not much 
difference between the weighted and unweighted analysis, which 
shows that the overall photometric quality (from observers HUZ 
and CTE) is quite homogeneous and there is not much effect from 
outliers. 

Fourier analysis of the middle season (160 days) showed a low 
frequency artefact (0.0044 c/d, a period of 222 days), a sampling 
artefact (0.996 c/d), a scatter of power between 2.7 and 2.99 c/d and 
again the 9.678 c/d peak. As with using the whole dataset, phasing 
the data with the doubled 0.68 d period produces a double humped 
light curve with slightly different amplitudes for each hump. 

Fourier analysis of the the last 30 days of data, which is the 
best subset, produced the results shown in the last three columns of 
Table[3] The 2.9 and 9.7 c/d frequencies show up in every analysis. 
The phase diagram at 2.9 (or 1.5) c/d shows a non-sinusoidal shape 
with sharp changes which is very suggestive of geometric varia- 
tions. The 1.52 c/d frequency may be a pure subharmonic within the 
uncertainties of the analysis and the 9.678 c/d peak is very strong. 

Rotational modulation can cause double humps if, for exam- 
ple, the modulation is caused by fixed star spots. So, if rotation is 
the cause, the formal frequency associated with the 2.9 c/d peak 
from the Fourier spectrum needs to be halved to obtain the true 
frequency. Also the 2.9 c/d folded light curve seems to be a bit 
different in shape when the whole dataset or different subsets are 
used. If this is due to rotation then shifts are possible if there is dif- 
ferential rotation on the B-star and/or spots change their latitudes. 
The 9.678 c/d frequency is very stable and very coherent (perhaps 
suggestive of some kind of geometric changes) and it is only 6 to 
6.5 times longer than the X-ray pulse period of 25 minutes. 




6 8 10 
Frequency (c/d) 

Figure 8. Periodograms obtained by using all of the HUZ and CTE data, (a) 
The original periodogram with the observing window shown in the inset, 
(b) Prewhitened periodogram showing the peak at 2.93 c/d. (c) Prewhitened 
periodogram showing the peak at 9.68 c/d. 



3.4 Radial Velocity Measurements 

Two approaches were taken to radial velocity determi nation: a cross 
correl ation with a model spectrum as computed by iMunari et al.l 
d2005l) and direct line profile fitting for the cores of the four domi- 
nant absorption lines. 

For the cross correlation approach, a high resolution (0.05 
A/pix) template was used after being degraded to the resolution 
of the data. The template spectrum (shown at the bott om of F ig.[4ll 
was fr om the library of model spectra as computed bv lMunari et al.l 
fcOOSl) , with the following parameters: T eff = 24000K (see Ta- 
ble Hi, log g = 3.0 (cgs units), [M/H]=-0.5, v rot = km 
s , microturbulent velocity = 2 km s~ , no [a/Fe] enhancement, 
old opacity distribu tion function (ODF) and no overshooting (see 
IMunari et al] d2005h for details). 

The observed spectra were cross-correlated with the tem- 
plate to compute relative radial velocities using the fxcor task 
of IRAF. Bary centric corrections were applied using the IRAF task 
rvcorrect. For direct line profile fitting, a sum of Gaussian and 
Lorentzian functions was first tried but ultimately a simple Gaus- 
sian fit of the line cores was used due to the limited S/N of the 
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Figure 6. Night by night presentation of season 2 V band data. Abscissas are HJD— 2450000, ordinates are instrumental V magnitude. Observer symbols are 
as in Fig. [2] in this case HUZ is the only observer. Some nights where only a few data points were taken are not shown. Note the presence of short period 
fluctuations (~0. Id) on every night. 
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Figure 7. Night by night presentation of season 3 V band data. Abscissas are HJD— 2450000, ordinates are instrumental V magnitude. Observer symbols are 
as in Fig. [2] Some nights where only a few data points were taken are not shown. 
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Table 3. Frequencies found from the iterative Fourier analysis of the HUZ and CTE data, listed in order of significance. Frequencies that were found in the 
third season data subset (the best one) plus at least one of the other two analyses are highlighted in bold face. 
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Figure 9. Photometric V band data from CTE folded two frequencies, (a) 
Folded at 2.924 c/d. (b) Folded at 1.462 c/d. The lower frequency appears 
to be composed of two peaks at slightly different amplitudes. 



Folded at period = 0.342 days 




Figure 10. Radial velocity data folded with a 0.342d period. The dashed 
line shows the best fitting sinusoid with a 0.342d period. 



data. Of the four dominant lines in the spectra, the weakest line (Si 
II at A4430) was dropped for all spectra and the remaining three 
lines were fitted. For the noisier spectra (due to clouds) the He I 
A4387 line could not be fitted, leaving only two fitted lines for 
those spectra. The resulting Doppler shifts were converted to helio- 
centric velocities using rvcorrect and the average and standard 
deviation of the line velocities were computed for each spectrum to 
give an uncertainty measure. The results (average line velocities) 
were comparable to the results obtained from cross-correlation. We 
considered the cross-correlation radial velocity values to be su- 
perior to the values obtained from individual line-fitting but the 
cross-correlation method gave no directly usable measures of error 
since the centreld peak-finding algorithm was used in IRAK 
The centreld algorithm was determined to be better for finding 
the maximum of the cross-correlation function than function fit- 
ting methods because the cross-correlation function was frequently 
asymmetrical. So we used the standard deviation (SD) of the line- 
fits to quantify the error in the cross-correlation radial velocity val- 
ues. The results, excluding the three noisy spectra are plotted in 
Fig. [TO] The radial velocity values for the individual lines were av- 
eraged here but will be useful individually for a larger data set be- 
cause each line potentially gives information about different layers 
in the Be star. 

Large radial velocity changes were seen in the course of a sin- 
gle night so the data were folded at the periods, and multiples of 
those periods, found in the period analysis of the photometry data. 
The radial velocity dataset is too small to permit any more exten- 
sive period searches. After the radial velocities computed from the 
three obviously noisy spectra were discarded, the best results were 
found at the prominent period of 0.342d. A sine curve was formally 
fit to the data folded at that period using a nonlinear least-squares 
Marquardt-Levenberg algorithm. The resulting fit, along with the 
measured radial velocity data are shown in Fig. QTJJ For the sine 
curve fit, Xred = 1-07 (5 degrees of freedom) where X^ e d i s tne 
reduced \ 2 value obtained by weighting the terms in the residual 
sum of squares with the SD values from the line-fits. The systemic 
velocity from the fit was 7 = —37.0 ± 4.3 km s _1 and the ra- 
dial velocity semi-amplitude was K\ = 21.8 ± 3.6 km s" 1 . We 
should emphasise that this conclusion about the radial velocity pe- 
riod is uncertain and subsequent data may lead to support for a 
period other than the ones tried here. 

The value of 7 = —37.0 ± 4.3 km s -1 closely matches th e 
radial velocities for othe r stars in NGC 663 dLiu et al.ll989l . ll99ll) . 
iMermilliod etail J2008l) report a cluster mean radial velocity of 
— 33.09±0.34 km s _1 with a dispersion of 1.74 km s" 1 . Our mea- 
surements therefore provide further confirmation that the HMXB 
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RX J0146.9+6121/LS I+61°235 is a member of the open cluster 
NGC 663 if we assume that the velocity kick from the supernova 
that formed the neutron star was non-existent or small. This is likely 
if RX J0146.9+6121/LS 1+61° 235 belong s to the class of BeXs 
with long orbital periods jPfahl et alj|2002l) as we suspect for this 
system (see the Discussion). One possible mechanism for such a 
small or non-existent supernova velocity kick is an electron-capture 
collapse of a degenerate O-Ne-Mg core, of a lower mass progen- 
itor, t hat can only happen in binary systems iPodsiadlowski et all 
120041) , 



4 DISCUSSION 

LS 1+61° 235 (V831 Cas ) was identified as t he optical counterpart 
of RX J0146.9+6121 by iMotch et all dl99ll) . They first identified 
the X-ray source from ROSAT data and subsequently observed Ha 
emission in the optical counterpart to confirm the BeX classifica- 
tion of the source. RX J0146.9+6121 is a persistent low-luminosity 
BeX with low X-ray variability (L m „/L mi „ ^ 10) an d it contains an 
X-ray pulsar with the longest pulse period known (Reig & Roche 

dMere 



Table 4. The astro-physical parameters of LS I+61°235 <Reig et al.ll997l) . 



1999) . The pulsar period is variable around 1400s (Mereghetti et al. 

2000) with variability presumably as the result of the increased 
mas s transfer a s sociat ed with occasional X-ray bursts. 

ICoe et al] d 19931) placed the distance at 2.2 kpc based on 
a reddening of E(B — V) = 0.92 deriv ed from uvby photom- 
etry (per the method of ICrawfordl dl978h ). This distance deter- 
mination is consistent with LS 1+61° 235 b eing a member of the 
open cluster NGC 663 l lPigulski et al]200ll) and was confirmed by 
iReig et al.l fl997h who further determ i ned th e astrophysical param- 
eters listed in Table l4l IPigulski et all d200lh give the age of NGC 
663 to be 20-25 Myr. The age of NGC 663 has been difficult 
to determine because the cluster contains the highest known Be 
star abundance in the Galaxy, and Be stars te nd to occupy anoma- 
lous positions in colour-m agnitude diagrams dFabregat & Torrei6r] 
120001) . iTapia et alj dl99 lb . for example, give the lowest determi- 
nation of 9 Myr. Usin g Stromgren uvby co l our in dicies for non- 
Be stars in the cluster, IFabregat & TorreiorJ dioooh deduce an age 
of 23 Myr in agreement wit h the determination by | Pigulski et al.l 
d200ll) . ICoe et all d 19931) and lReig et all d 19971. l2000l) obtained IR 
JHK photometry that showed long term variations of approximately 
0.3 mag for J (J~9.9), 0.3 mag for H (H~9.6) and up to 0.5 mag 
for K (K~9.4). 

The Be decretion disk of LS I+61°235 has been intensively 
studied through the investigation of Violet/Red (V/R) variations in 
the Ha emission. The Ha, Hel A6678 and the Paschen lines Pal 1, 
A886 3, and Pal2, A8750 all show variable emission dReig et all 
l2000h . Furthermore, the Ha emission shows a quasi-cyclic varia- 
tion of 1240 ± 30 days. The V/R variation of the Ha line corre- 
lates with long-term change s in Ha EW and IR (JHK) intensity. As 
shown bv lReig et al.l d2000l) . the V/R and IR variat ions can be ex- 
plained by the global one-arm ed oscillation model dOkazaki|[l99ll , 
ll997l : |Papaloizou et alj[l992T) consisting of global m = 1 oscilla- 
tions of the cool Be disk in which an enhanced density perturbation 
develops on one side of the disk and slowly precesses, in this case 
in the prograde direction, with one precession cycle being equal to 
the V/R cycle. The density perturbation is confined to a few stellar 
radii in the disk and the precessio n period is expected to be rel - 
atively insensitive to the disk size dSavoniie & Heemskerkiri993l) . 
Since the variations are azimuthal, no variation in IR colour is ex- 
pected and this is what has been observed. Radial variation in the 
disk would cause IR colour variation. We will relate the motion of 



Spectral type 


B 1 V 


E(B - V) 


0.93 ± 0.02 


Terr 


24000 ± 1500 K 


Radius 


7±1R 


Mass 


11 ± 2M 


My 


-3.1 ±0.5 


Distance 


2.3 ± 0.5 kpc 


BC 


-2.4 ± 0.2 


M bo , 


-5.5 ±0.5 (1.2 x 10 4 L Q ) 


log g 


3.9±0.21og(cm/s 2 ) 


v m i sin i 


200 ±30 km s" 1 



the one-armed density perturbation with the X-ray behaviour of RX 
J0146.9+6121 below. 

The strongest results arising from our data for LS 
I+61°235/RX J0146.9+6121 are three frequencies in the V band 
light curve: 2.9 c/d, 1.5 c/d and 9.7 c/d. That is, the periods 0.34d, 
0.67d and O.lOd are present. Such short periods have been argued 
for other Be stars to be due to rotation or non-radial pulsation 
dBalomJl990l) . The 0.67d period is almost surely a subharmonic of 
the 0.34d period and, if the light variation is caused by a geomet- 
rical effect such as ellipsoidal variation (the Z-shape of the folded 
light curve implies a more complex geometry than simple ellip- 
soidal variation), then the spin period of the Be star is 0.68d. This 
conclusion can be checked using the astrophysical par ameters for 
LS I+61°235 previously published bv lReig et all d 19971) (see Table 
|4j. Using the value of 2.9238 c/d (0.34d period) as determined from 
the analysis of all three observing seasons and a Be star radius of 
7±1 Rq we find a surface velocity of 1036 ±148 km s _1 . With a 
rotation frequency of half that, a surface velocity of 518 ±74 km 
s _1 results. The break-up velocity is 



GNU 
R* 



(1) 



Using NU = 11 Mq and R, = 7 Rq, the critical velocity is 
D cr j t = 547 km s _1 . So the rotation period cannot be 0.34d or 
the star would break-up. A rotation period of 0.68d is near-critical, 
a condition that is b elieved to hold for Be s tars (e.g. from statistical 
studies like that of IChauville et al .1 d200ll) . from the direct inter- 
fer ometric measurement of the oblat eness of the Be star Achernar 
bv lDomiciano de Souza et al] d2003h and, from theoretic al consid- 
eration s on how the decretion disk is formed like that of iLee et"al] 
dl99ll) ). So it is possible that the rotation period may be 0.68d with 
the 0.34d period representing when opposite sides of the star are 
presented to us. 

With a rotation period of 0.67d and v ro t — 518 km s _1 we can 
use the previously determined value of v m t sini = 200 km s _1 to 
deduce that the inclination is i — 23°. Taking into account the error 
reported for R* and v m t sin i we have i — 23 _ 6 degrees. 

Next we investigate the possibility that the 0.34d period is the 
binary orbital period since we have some evidence that the radial 
velocity period is 0.34d. The mass function, /(M), with masses 
Mi, the mass of the Be star and M%, the mass of the compact com- 
panion, given in solar masses, the orbital period P or (, in days and 
the semi-amplitude of the radial velocity variation, K\, given in km 



f(M) = 



(Mi ± M 2 



= (1.0361 xHr 7 )(l-e 2 )3A-iP orb (2) 
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where e is the eccentricity of the orbit (Hilditch 2001). Taking e = 
0, Mi = 11 M , A'i = 22 km s" 1 and i = 23° (assuming the Be 
spin axis and the orbital axis to be aligned) in Eq.[2]gives the mass 
of the compact companion as M 2 = 0.45 Mq for P or (, = 0.34d. 
The semi-major axis of the orbit is given by Kepler's law as 



^orb 

4-7T 2 



G(Mi + M 2 



(3) 



Using M 2 = 0.45 M Q and P orb = 0.34d gives a = 4.6 R Q . Since 
the radius of the Be star is 7±1 M , the proposed orbit is clearly 
too small, so 0.34d cannot be the orbital period. 

The mass and luminosity of LS 1+61° 235 as given in Table [4] 
put it in the middl e of the j3 Cep instability strip on the Hertzspru ng- 
Russell diagram JPamvatnvkblll999l;IStankoy & Handle j|2005l) . In 
a recent catalogue dStankov & Handlenl2005l) . f3 Cep stars are de- 
fined to be massive nonsupergiant variable stars with spectral type 
O or B whose light, radial velocity and/or line profile variations are 
caused by low-order pressure and gravity mode pulsations. For low 
order radial modes the light curve and radial velocity will have the 
same period. This is what appears here if we take the radial velocity 
period to be 0.34d. In that case the pulsation constant, Q, given by 



Q 



3M 
4TvR :i 



(4) 



where P is the period and M and R are the mass and radius of the 
star in solar units, is 0.03±0.01d, if the astrophysical parameters 
from Ta ble B]are used. The moda l Q of known /3 Cep variables is 
0.033d dStankov & Handlej|2005h which corresponds to the radial 
fundamental mode. 

If, between the two periods 0.34d and 0.67d, we take 0.34d 
as the pulsation period, we can not interpret the 0.67d period as a 
harmonic due to a non-sinusoidal pulsation in time. Thus the inter- 
pretation of the 0.67d period as the Be star spin period with a spin 
axis inclination of 23° may still hold. The combination of 0.34d be- 
ing a harmonic (or close to a harmonic) of the 0.67d period plus a 
pulsation period of 0.34d would make the 0.34d period the promi- 
nent one, as observed, since two processes would contributing to 
the light variation at that period. A simpler interpretation is that the 
0.34 radial velocity period is spurious and the 0.34d photometric 
period is simply a harmonic of the 0.67d period due to the non- 
sinusoidal shape of the 0.67d period. Further radial velocity data 
can resolve this ambiguity. If the O.lOd period is due to pulsation 
and 0.34d is the fundamental mode, then it must be a higher order 
p mode. 

Finally, we speculate on the origin of the "super" X-ray out- 
bursts spaced at approximately 450d. The prograde period of the 
one-armed density enhancement in the Be disk h as been measure d 
to be approximately 1240d with some variability dReig et aljioooh . 
If we assume that the compact star also orbits in a prograde direc- 
tion (relative to the Be star spin) then we might expect a super out- 
burst to occur when the arm and the compact star line up. Given the 
arm period of 1240d and an outburst (alignment) period of 450d, 
this means that the compact star's orbital period is 330d in agree- 
ment wi th an expectation of abou t 300d from the Corbet spin-orbit 
relation dCorbetll986l : lReig||2007h . This conclusion is roughly con- 
sistent with the 184d period found in the RXTE ASM data which 
would be explained by the occurrence of X-ray bursts when the 
compact star passed through the plane of the decretion disk. The 
cause of the apparent 90d intervals between the minor outbursts, if 
these intervals are significant, is, however, not clear in this interpre- 
tation since the minor outbursts would have to occur four times per 



orbit. Perhaps, if the orbit is sufficiently inclined to the Be disk, the 
mass flow onto the neutron star will switch poles 90° in the orbit 
from the disk crossing nodes and this switch, through some as yet 
undetermined mechanism, could lead to a minor X-ray outburst. 

An orbital period of 330d matches expectations from the Cor- 
bet spin-orbit relation but is it consistent with the X-ray luminosity 
of RXJ0146. 9+6121 an d typical mass loss r ates for Be stars? Based 
on a distance of 2.5 kpc. lHaberl et al. I dl998l) determined that the X- 
ray luminosity, L x varies between 2 x 10 34 and 2 x 10' !5 erg s _1 . 
The X-ray lumiosity is related to the mass accretion rate, M 2 onto 
the neutron star through the conversion of gravitational potential 
energy to luminosity by 



M 2 = 



L X R 2 
GM 2 



(5) 



where G is the gravitational constant and where we may take 
R2 — 10 6 cm and M2 = 1.4 M to be the radius and mass of 
the neutron star. So M 2 is between 2 x 10~ 1 2 and 2 x 10" 11 M 
yr -1 . Following lFrank. King & Rainei d2002l) we assume that mass 
leaves the Be star at escape velocity and will be captured within a 
cylindrical region that has an axis through the neutron star in the 
direction of the relative wind velocity. Defining v w to be the wind 
velocity at a large distance from the Be star and v IC \ as the wind 
velocity relative to the neutron star, we have v w = v2« cr it (the 
escape velocity, see Eq.[T]for the definition of w C rit) and v IC \ = v w 
since the Keplerian velocity about the Be star at the neutron star 
distance, a, is small relative to v w . So the cylindrical capture axis 
direction would be essentially towards the Be star. The radius of the 
cylindrical capture region, r- acc , is where the kinetic energy of the 
wind equals the gravitational potential of the neutron star dBondil 



2gm 2 



(6) 



The ratio of the emitted wind, —Mi, to the captured wind, M 2 , 
mass flow rates is therefore 



(7) 



-Mi 47ra 2 



where the binary separation, a, is given by Eq.[3] If v Ic i = \/2v CT it 
then 



_M 

^M; 



1 /MaA 2 (Ri 
1 4 VMj I a 



(8) 



For Mi = 11 M , M 2 = 1.4 M and an orbital period of 330d, 
a = 480 Rq, and using Ri = 7 R in Eq. [8] gives a ratio of 
emitted to captured wind of 10~ 6 . So the inferred Be wind rate, 
-M, would be between 2 x 10" 6 and 2 x 10" 5 M yr" 1 . This 
is consistent with Be mass outflow rates of ~10~ 6 M yr -1 seen 
observationally and exp ected on theoretical grounds dLamers et al.l 
l200d . lVink et al .1 1 2000). However, we should also note that some 
observational and theoretical studies have come to the conclusion 
that the Be mass lo s s rate is much lowe r at 10 -9 to 10 -8 M yr -1 
dWaters et al.lll988l ; lHaberl et al .11 1998b . 



5 CONCLUSION 

The V band light curve of LS 1+61° 235 shows three strong period- 
icities at 0.34d, 0.67d and 0. lOd. In our interpretation, the 0.34d 
period may be from pulsation in the radial fundamental mode, 
the 0.68d period is the Be star spin period and the 0. lOd period 
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is a higher order p mode pulsation. The spin period, when com- 
bined with a previous determination of the astrophysical param- 
eters of the Be star, gives an inclination for the Be spin axis of 
i = 2 3 degrees. Comparing the latest catalogue of /3 Cep vari- 
ables JStankov & Handlenl2005t) to the latest catalogue of Galactic 
HMXBs dLiu et alj|2006l) reveals that LS 1+61° 235 would be the 
first confirmed /3 Cep member of a HMXB if subsequent observa- 
tions uphold these conclusions. The systemic radial velocity con- 
firms that LS 1+61° 235 is a member of the cluster NGC 663. If we 
interpret that the X-ray "super" outbursts happen when the compact 
star lines up with the one-armed Be decretion disk enhancement, 
then the orbital period is approximately 330d. This orbital period 
agrees with the value exp ected from the Corbet spin-orbit relation 
< ICorbetlll98d : lReidl20fA 

The variability of the Ic light curve looks to be similar to the 
variability of the V band light curve and there may be some vari- 
ability in B — V colour. More observations are needed to investi- 
gate these possible colour variations which should exist if pulsation 
is present. 
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